Jamming Percolation and Glass Transitions in Lattice Models 
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A new class of lattice gas models with trivial interactions but constrained dynamics are introduced. 
These are proven to exhibit a dynamical glass transition: above a critical density, p c , ergodicity is 
broken due to the appearance of an infinite spanning cluster of jammed particles. The fraction 
of jammed particles is discontinuous at the transition, while in the unjammed phase dynamical 
correlation lengths and timescales diverge as exp[C(p c — Dynamic correlations display two- 

step relaxation similar to glass-formers and jamming systems. 



In the majority of liquids dramatic slowing down oc- 
curs upon super-cooling. In a rather small temperature 
window, typically from the melting temperature T m to 
about 2T m /3, the viscosity increases by 14 orders of mag- 
nitude and the relaxation becomes complicated: non- 
exponential and spatially heterogeneous Similar fea- 
tures are observed in soft materials, such as colloidal sus- 
pensions and more generally in non-thermal "jamming" 
systems 0- Despite a great deal of effort, these remark- 
able phenomena, associated with "glass transitions" are 
still far from understood. Even the most basic issues re- 
main open: is the rapid slowing down due to proximity 
to a phase transition? Is this putative glass transition 
static or purely dynamic 0]? Experimental results make 
it clear, however, that if an ideal glass transition does 
exist it should have some peculiar features: the density 
auto-correlation function C it) should exhibit a lengthen- 
ing plateau that, at the transition, extends out to infinite 
times. Thus, the Edwards-Anderson order parameter, 
defined as the infinite time limit of C(t), will be discon- 
tinuous at the transition. But this discontinuity should 
be accompanied by a critical divergence of the relaxation 
time. And, contrary to usual critical slowing down, the 
relaxation time should diverge exponentially, as in the 
Vogel-Fulcher law Long-range spatial correlations, if 
they exist at all, must be very subtle. These unusual 
properties present major theoretical challenges: whether 
or not there is a true transition there is no "standard" 
framework to start from. There are promising results for 
models on Bethe lattices j4( and for some with long-range 
interactions |5| . But the quest for models with short range 
interactions and no quenched disorder that are simple 
enough to be analyzed and can be shown to have a glass 
transition, namely a transition with the basic properties 
discussed above, is still open, in spite of much effort. 

In thispaper we introduce the first examples of such 
models [fj, Ijj. These are kinetically constrained mod- 
els (KCMs): stochastic lattice gases with no static in- 



teractions, except hard core, but constrained dynamics 
. The elementary moves are particle jumps for conser- 
vative models and birth/death moves for non conserva- 
tive models. Whether a move can occur depends on the 
nearby configuration and is non-zero only if some local 
constraints are satisfied. These kinetic constraints can 
radically change the dynamical behavior and typically 
induce glassy phenomenology |il 1^. Hfl| . For some KCMs 
the dynamics becomes so slow at high density or low tem- 
perature, that they have been conjectured to undergo a 
true glass transition. The simplest examples are Kob- 
Andersen models on a square lattice (SKA) where 
particles can move if and only if they have no more than 
two nearest neighbors both before and after the move. 
Although the analogous model on a Bethe lattice 0, 0] 
does have a jamming transition, we have shown previ- 
ously in 0| that the SKA and a broad class of gener- 
alizations on hypercubic lattices cannot have ergodicity 
breaking transitions: in any finite dimension the relax- 
ation time diverges — in many cases in a super- Arrhcnius 
way — only at the close packing density (p = 1). 

But this is not the only possible behavior. We here in- 
troduce a new class of KCMs that do exhibit a jamming 
transition at a non-trivial critical density, p c on finite di- 
mensional lattices. 

For simplicity we focus on one of the simplest: a square 
lattice model without particle conservation; vacancies can 
loosely be thought of as "free volume" which need only 
be conserved on average. At the end, we discuss gen- 
eralizations to both higher dimensional and conservative 
models. The stochastic dynamics is as follows: An oc- 
cupation variable at site x cannot change if x is blocked 
along either of the diagonal directions, as defined below. 
Unblocked sites change from occupied to empty and from 
empty to occupied with rates (1 — p) and p, respectively. 
Thus detailed balance is satisfied with the trivial product 
measure with density p. The blocking is determined by 
the eight fourth-nearest-neighbor sites of x. Denote pairs 
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FIG. 1: a): Sites connected by arrows form one of the graphs 
on which directed percolation can form frozen clusters: e.g. 
the set of occupied sites (big dots) shown. b)Site X and its 
NE,NW,SW,SE pairs of neighbors, c) Portion of an empty 
octagonal annulus. The interior of the annulus, as any fi- 
nite region surrounded by a double frame of vacancies, can 
be eaten away. Whether the vacant region can expand is de- 
termined by the three key sites indicated by question marks. 
If one of these belongs to an occupied DP path in the NE 
direction which is anchored on two perpendicular DP paths, 
as shown, it blocks the expansion. Inset: the full octagon. 
A necessary condition for the octagon not to be expandable 
of one step in the NW direction is that a DP path spans the 
dotted rectangle. 

of these the north-east (NE), north-west (NW), south- 
east (SE), and south-west (SW) pairs as in Fig. Q]b). 
Site x is blocked if either at least one of the NE sites 
and at least one of the SW sites is occupied, or at least 
one of the SE sites and at least one of the NW sites is 
occupied. Blocking can thus be along either the NW-SE 
or the NE-SW diagonals. As the distance to the blocking 
neighbors resembles a knight's move in chess, we call this 
the "knights model". 

If a site cannot be unblocked even by first emptying with 
allowed moves an arbitrarily large number of other sites, 
we call the site frozen. Any finite cluster of particles 
cannot be frozen: one can always unblock all sites by 
emptying from the perimeter in (see Fig. ^i). A crucial 
question is whether an infinite spanning cluster of frozen 
sites exists in infinite systems. We will call this problem 
jamming percolation; it is akin to bootstrap percolation 

We will show the following results (which can be proved 
0]): (i) w ith blocked or periodic boundary conditions 
on L x L squares, there exist configurations with system- 
spanning clusters of frozen sites; (ii) on infinite lattices, 
below a critical density, p c , there are no infinite frozen 
clusters; while (iii) above p c , there is an infinite cluster 
of frozen particles that occupies a non-zero fraction, (p^, 
of the area; (iv) p c coincides with the critical density 



for directed site percolation (DP) on a square lattice; 
(v) 4>oo is discontinuous at p c \ (vi) below p c , there is a 
crossover length S(p): squares of size L << S are very 
likely to have a frozen cluster, while for L >> S, the 
probability of a frozen cluster falls-off exponentially; (vii) 
as p increases to p c , 5 diverges exponentially rapidly, as 
log 5 ~ (p c — p)~ M with fi = 0.64 related to DP exponents; 
(viii) the relaxation time diverges as 5 or faster. Thus 
even though the critical density is the same as for DP, 
the behavior is completely different. 

We first show that both an unfrozen and a frozen phase 
exist. At sufficiently low densities, the occupied sites 
will not percolate via connections up to fourth neighbors: 
the resulting finite clusters can always be unblocked from 
their perimeters. Thus at low densities, frozen clusters 
will not occur. In contrast, at high densities, spanning 
frozen clusters occur. Consider site directed percolation 
with directed links that connect a site to its two NE 
(fourth) neighbors as in Fig. [I] a). Infinite directed paths 
of occupied sites exist for p > p^ p = 0.705 (the criti- 
cal threshold for conventional site DP on a square lattice 
[l5]|). These clusters of sites are frozen since each has at 
least one occupied fourth-neighbour in both NE and SW 
directions. Thus, p c < p® p . 

For the above argument it was sufficient to use block- 
ing along just one of the two diagonal directions. But, 
due to blocking in the perpendicular diagonal direction, 
typical frozen configurations do not resemble DP clus- 
ters: they consist of short DP paths that terminate at 
each end in a T-junction with a DP path in the perpen- 
dicular direction. Thus large regions can be frozen even 
if they are not spanned by DP clusters. 

An explicit construction is instructive. Consider a 
structure built of DP paths of length, s, intersecting at 
T-junctions as in FigJ^Ja). This structure does not con- 
tain any long DP cluster, yet it is frozen. And, crucially, 
there exists a similar frozen cluster as long as each DP 
path remains inside a nearby rectangular region of size 
sxs/6 (see Fig. [21 a) - Therefore the probability of the 
system being frozen is bounded from below by the prob- 
ability that all these rectangles are spanned. This will be 
substantial if the DP spanning probability of each such 
rectangle is very high. 

The crucial needed property follows from the 
anisotropy of critical DP clusters: a cluster of length s 
typically has transverse dimensions of order wit h C , 
(often called l/z) the anisotropy exponent, £ = 0.63 [15j. 
Therefore an s x bs rectangle with s in the parallel direc- 
tion can be cut into slices of width such that for each 
slice the probability of DP spanning paths is substantial 
for p > p® p . Therefore, even at p® p , the probability 
rt>(s) of not having any DP crossing in a large s xbs rect- 
angle is rb(s) < 0[exp(— 2Cbs 1 ~ < >)] (with C a constant). 

What happens just below p,P p ? The above argument 
will still hold for rectangles that are are of order the 
DP parallel correlation length, £m. This, together with 



the previously explained construction of a frozen struc- 
ture, implies that an L x L square is likely to have 
a frozen structure built of DP clusters of length £|| if 
ri/^)L 2 /q « 1. For L < S< ~ exp(C<£}- c ) - 
a lower bound on S — squares will thus contain frozen 
clusters with high probability in spite of the rarity of DP 
clusters with length larger than £|| <C S . 

We now need to show that below p® p , sufficiently large 
squares are unlikely to contain frozen clusters. Again, 
this can be done by construction — now of unfrozen re- 
gions. Consider an infinite system within which is an 
octagonal annulus of radius (center to NW) I that is com- 
pletely empty. Whether or not this empty region can be 
expanded depends crucially on three key sites in each di- 
agonal corner, say, the NW corner, as shown in Fig. ^ 
c). If all three key sites are empty or emptyable — i.e. 
unfrozen — the empty annulus can be expanded along 
its two adjacent (NNW and WNW) sides. A necessary 
condition to have a key site frozen is that it belongs to a 
DP cluster in the NE direction that is anchored at both 
ends, as in Fig. He). Furthermore, in order for this an- 
chorage to occur, it is necessary that the NE path spans 
the rectangular dotted region of size £ x b£ (with b a con- 
stant) in Fig 1 c). For p < pf p DP clusters with length 
£ much larger than the DP correlation length, £|| , are ex- 
ponentially rare and the rectangle spanning probability 
is - exp(-£/£||) 0. Thus, if p < pf p , the probabil- 
ity that the annulus can be expanded out to infinit y by 
successive expansions is high for £ » £|| ||, Ull Il3|. 
Since the infinite system will contain a non-zero density 
of these empty regions which can be expanded to unblock 
the whole system, we conclude that p c = p® p . 

Estimating how rare the unblocking regions are near 
Pc yelds an upper bound on the crossover length S. 
Starting with a small empty nucleus the — small — prob- 
ability, (5, that it can be expanded out to size ~ £|| (and 
hence readily to infinity) is the product of many small 
terms. This is dominated by £ ~ £||: from f(£||)i we 
obtain 5 > exp(— 2C>£|j ). Since in an L x L square 
there are (L/^) 2 roughly independent positions for such 
empty nuclei, some are likely to occur if 5L 2 /^ 2 is not 

small. Thus we find that S < S> ~ £|| exp(C>£|j~ 1 '). 
We have found upper and lower bounds for the crossover 
length, S, of similar form, hence 

logS~£^~ (p c -p)->* with (j, = (1 - C)v\\ = 0.64 

with ^|| = 1.73 the correlation length exponent for DP 

We now discuss another peculiarity of this transition: 
the nature of the frozen clusters implies that the density, 
4>txi °f the infinite frozen cluster jumps at p c . To analyze 
the probability that a site is frozen, consider an occu- 
pied site, e.g. the origin, which belongs to a DP cluster 
that extends to a (small) distance £q/2 in both the NE 




FIG. 2: a) Frozen structure made of DP paths of length s 
(represented by straight continuous lines) intersecting at T- 
junctions. The A path anchors one end of the B and C paths, 
while its ends are anchored by the D and E paths. This an- 
choring is retained if path A is displaced until the nearby 
dotted line, even if the B and E paths are similarly displaced 
no further than their corresponding dotted lines. Thus the 
structure is frozen if all the rectangles — shown and unshown 
— formed by the solid and dashed lines are spanned length- 
wise by DP paths b) Two sequences of intersecting rectangles, 
IZi, that when spanned length- wise by directed paths, make 
the central site, O, frozen. For clarity, the rectangles in one 
direction are shown with dashed lines. 



and SW directions: this occurs with some probability pq. 
Now focus on two infinite sequences of rectangles IZi of 
increasing size £j x £i/12 with £\ — £o, £i — 2£;_2 and 
intersecting as in Fig. [2J If each of these rectangles is 
spanned lengthwise by a DP cluster, the origin is frozen. 
The probabilities of perpendicular DP paths are posi- 
tively correlated. Thus <\>oo > po Iljl — r i/i2(^i)] 2 - Since 
£|| is infinite at p c , as shown above, rt>(£) decays expo- 
nentially to zero as exp[— Cb£ Y ~ z \. Therefore, the infinite 
product is non-zero, giving a strictly positive bound on 
^oo(pc) for any £q. Thus, in contrast to DP or conven- 
tional percolation, the infinite frozen cluster of jamming 
percolation is "compact" — i.e. with dimension d — at 
the transition. 

To supplement our predictions, we have studied L x L 
systems numerically. The distribution of the densities, 
(j)L, of the frozen clusters shows two peaks with weak 
size dependence as found at conventional first order tran- 
sitions. This is consistent with the predicted discontinu- 
ous behavior. The probability that there exists a frozen 
cluster is substantial for p fifteen percent below p c even 
in our largest systems, (L = 1600): it is thus hard to 
study the asymptotic critical behavior (see [16j for an 
analogous problem in the context of bootstrap percola- 
tion) . But in a slightly different model one can get closer 
to the transition [lj|: these data are consistent with the 
predicted InS ~ (p c — p)~^ with \l = 0.64, but the small 
range of In L available makes the uncertainties in p large. 

Our results on jamming percolation have important 
implications for the equilibrium dynamics of the knights 



model. For p < p c the fact that infinite frozen clus- 
ters do not exist i mply that the system is ergodic [proof 
can be done as in (11] 2.5)] but dynamical correlations, 
such as the density auto-correlation function, C(t), dis- 
play two step temporal relaxation with a long plateau 
followed eventually by a relaxation (numerical results 
will be reported elsewhere [13 )• The relaxation time di- 
verges exponentially near p c , at least as S z , with z > 1: 
This is reminiscent of the Vogel-Fulcher law found exper- 
imentally near glass transitions. Above p c , the plateau 
stretches to infinite times with the Edwards-Anderson 
order parameter, q = limt^oo C(t), discontinuous at the 
transition. This follows from our results, since q is related 
to the density of frozen sites. 

We have seen that, even though the critical densities 
are the same, the properties of jamming percolation are 
strikingly different from the power law behavior of di- 
rected percolation. Most of the physics is controlled 
by relatively short DP clusters joined together at T- 
junctions. The only role of long DP clusters is to prevent 
very rare large unfrozen regions from unblocking their 
surroundings. There is substantial universality in the pri- 
mary features of the jamming percolation. This extends 
even to the SKA model which has frozen configurations 
composed of double-width occupied bars that terminate 
in T-junctions with similar perpendicular bars, but there 
is no real transition because very long bars are unlikely. 
Yet the SKA's behavior as p f 1 is similar to the knights 
model as p f p c with 1/(1 — p) roughly replacing powers 
of £|| . Surprisingly, if the square lattice of the SKA is re- 
placed by a particular complicated four-fold coordinated 
lattice, p c < 1 and the behavior is similar to the knights 
model. Thus the local structure matters a lot as in real 
glasses. Note that cooperative models with a transition, 
in contrast to those without, display two-step relaxation 
from the dynamics within blocked regions: this is like 
beta-relaxation in glasses . 

Models with particle- conserving dynamics behave sur- 
prisingly similarly to those without: the nature of the 
transition (and in some cases the critical density) re- 
main the same because the slowing of the dynamics is 
dominated by the large clusters of the underlying jam- 
min g pe rcolation. Diffusive transport rides on top of this 
[hi. Il7l Hif. In three dimensions, two natural general- 
izations of our jamming percolation exist: one composed 
of DP clusters — which should slow down as a double 
exponential of (p — p c )~^ — and the other of directed 
sheet-like structures which will have exponential slowing 
down like we have found in 2D. The key ingredients are 
kinetic constraints that enable huge jammed clusters to 
form out of small objects without these becoming much 
more common or much larger. 

For the future, the connection between our results and 
the jamming transition found for continuum particle sys- 
tems needs exploring. With the hope of increased 
understanding of the rapid liquid to glass crossover ob- 



served experimentally, one should also analyze the effects 
of constraint-violating processes occurring with a very 
low rate. For both glasses and granular materials, study- 
ing the non-e quili brium effects caused by a quench or by 
driving forces |14j is merited even in the simplest models 
that exhibit a jamming transition. 

After the completion of this work a new version of a 
preprint Q appeared in which other models with a jam- 
ming transition are introduced and studied numerically. 
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